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1.1 Introduction 

The control of transport in magnetically confined plasmas is of major importance 
in the long way to achieve controlled thermonuclear fusion. Two major mechanisms 
have been proposed for such a turbulent transport: transport governed by the fluc- 
tuations of the magnetic field and transport governed by fluctuations of the electric 
field. There is presently a general consensus to consider, at low plasma pressure, that 
the latter mechanism agrees with experimental evidence [1] . In the area of transport 
of trace impurities, i.e. that are sufficiently diluted so as not to modify the electric 
field pattern, the E x B drift motion of test particles should be the exact transport 
model. The possibility of reducing and even suppressing chaos (combined with the 
empirically found states of improved confinement in tokamaks) suggest to investi- 
gate the possibility to devise a strategy of control of chaotic transport through some 
smart perturbations acting at the microscopic level of charged particle motions. 
Chaotic transport of particles advected by a turbulent electric field with a strong 
magnetic field is associated with Hamiltonian dynamical systems under the approxi- 
mation of the guiding center motion due to E x B drift velocity. For an appropriate 
choice of turbulent electric field, it has been shown that the resulting diffusive trans- 



1 



2 CHAPTER 1 . RED UCTION OF THE CHAOTIC TRANSPORT OF IMPURITIES IN TURB ULENT MAC 




X 



Figure 1.1: The interaction between the electric turbulent field E and the toroidal 
magnetic field B produces the E x B motion which causes a drift of particles from 
the center of the plasma toward the edge. The guiding center motion of charged 
test particles on a transversal plane to the toroidal direction can be described by a 
Hamiltonian formalism with the spatial coordinates (x, y) which represents canonical 
conjugate variables and with the electrostatic potential V, E = — VV, which plays 
the role of the Hamiltonian, as described by Eq. (11. ip . 

port is found to agree with the experimental counterpart [21 E] ■ 

Here we address the turbulent transport of particles governed by the interaction 
of the electric turbulent potential E generated by the plasma itself and the strong 
confining magnetic field B produced by external coils. The resulting E x B drift 
motion is perpendicular to the confining magnetic field and causes losses of particles 
and energy from the inner region to the border of the device, with a consequent 
decrease of plasma temperature (see Fig. 11. ip . The equations of this drift motion 
for charged test particle, in the guiding center approximation, have a Hamiltonian 
structure and are given by 

where x = (x, y) represents the spatial coordinates of the transversal section to the 
confining toroidal magnetic field, B the norm of the magnetic field B, c the velocity 
of light and V is the electric turbulent potential, that is E = — VV. We notice that 
the resulting dynamics is of Hamiltonian nature with a pair of canonically conju- 
gate variables (x, y) which consists of the position of the guiding center. Since the 
potential which plays the role of the Hamiltonian is time-dependent, it is expected 
that the dynamics is chaotic (with one and a half degrees of freedom). 
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The sensitivity of chaotic systems to small perturbations triggered a strong in- 
terdisciplinary effort to control chaos [U El El El El EEl EH] ■ After the seminal work 
on optimal control by Pontryagin [TS] , new and efficient methods were proposed for 
controlling chaotic systems by nudging targeted trajectories [131 ESI HI [T71 [TS]. 
However, for many body experiments like the magnetic confinement of a plasma or 
the control of turbulent flows, such methods are hopeless due to the high number of 
trajectories to deal with simultaneously. 

Here the control is performed with the addition of a control term / to the electric 
potential. It consists as a small and apt modification of the electric potential with 
relevant effects on the reduction of chaotic particle diffusion. More generally, con- 
trolling Hamiltonian systems here means to achieve the goal of suppressing chaos 
when it is harmful as well as to enhance chaos when it is useful, by driving the 
dynamical evolution of a system toward a "target behavior" by means of small and 
apt perturbations. For magnetic confinement devices for controlled thermonuclear 
fusion, the two possibilities of increasing chaos and suppressing it are of interest. In 
fact, in the core of a magnetic confinement device the aim is to have a stable and 
regular dynamics in order to induce as many fusion reactions as possible, that is 
a high rate of energy production. At the same time, in order to spread heat over 
a large area, which is fundamental both for a good recycling of energy and for the 
preservation of plasma facing components, at the edge of the confinement device the 
problem of draining highly energetic particles requires the capacity of increasing the 
degree of chaoticity of the dynamics. 

For reducing chaos, KAM theory gives a path to integrability and the possibility 
of controlling Hamiltonian systems by modifying it with an apt and small perturba- 
tion which preserves the Hamiltonian structure. In fact, if on one hand it is evident 
that perturbing an integrable system with a generic perturbation gives rise to a loss 
of stability and to the break-up of KAM tori, on the other hand the structure of the 
Hamiltonian system changes on the set of KAM tori with a continuity C°° [19] with 
respect to the amplitude of the perturbation. The idea is that a well selected per- 
turbation instead of producing chaos makes the system integrable or more regular, 
recovering structures as the KAM tori present in the integrable case. 

Starting from these ideas we have addressed the problem of control in Hamilto- 
nian systems. For a wide class of chaotic Hamiltonians, expressed as perturbations 
of integrable Hamiltonians, that is H + eV, the aim is to design a control term / 
such that the dynamics of the controlled Hamiltonian Hq + eV + / has more regular 
trajectories (e.g. on invariant tori) or less chaotic diffusion than the uncontrolled 
one. Obviously / = — eV is a solution since the resulting Hamiltonian is integrable. 
However, it is a useless solution since the control is of the same magnitude of the 
perturbation while, for energetic purposes, the desired control term should be small 
with respect to the perturbation eV. For example the control term should be of 
order e 2 or higher. In many physical situations, the control is only conceivable 
or interesting in some specific regions of phase space where the Hamiltonian can be 
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known and/or an implementation is possible. Moreover, it is desirable to control the 
transport properties without significantly altering the original setup of the system 
under investigation nor its overall chaotic structure. The possibility of recreating 
specific regular structures in phase space such as barriers to diffusion, can be used 
to bound the motion of particles without changing the phase space on both sides of 
the barrier. 

In addition, the construction of the control term has to be robust. In fact, one 
expects from a KAM theorem also a given robustness with respect to modifications 
of the well selected and regularizing perturbation. This fact still comes from the 
previous results that we mentioned above about the transition from integrability to 
chaos that Hamiltonian systems generically show. In other words, if an apt mod- 
ification of a chaotic Hamiltonian gives rise to a regularization of the dynamics, 
with for example the creation of KAM tori, one expects that also some approxima- 
tions of this apt modification are still able to perform a relevant regularization of 
the system. This is a fundamental requirement for any kind of control scheme in 
order to guarantee the stability of the results and in view of possible experimental 
implementations. 

1.2 Control term for Hamiltonian flows 

In what follows, we explain a method to compute control terms of order e 2 where e 
is the amplitude of the perturbation, which aims at restoring the stable structures 
which were present in the case where e — 0. We notice that such stable structures 
are present for Hamiltonian systems H + e'V where e' < e c < e (where e c is the 
threshold of break-up of the selected invariant torus), even though they are deformed 
(by an order e') . 

For a function H, let {H} be the linear operator such that 

{H}H' = {H,H'}, 

for any function H', where {• , •} is the Poisson bracket. The time-evolution of a 
function V following the flow of H is given by 




= {H}V, 



which is formally solved as 

V(t) = e*< H >y(0), 

if H is time independent. Let us now consider a given Hamiltonian H Q . The operator 
{H } is not invertible since a derivation has always a non-trivial kernel (which is 
the set of constants of motion). Hence we consider a pseudo-inverse of {Ho}. We 
define a linear operator T such that 

{Ho} 2 T = {H }, (1.2) 
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i.e. 

W, {H ,{H ,TV}} = {H ,V}. 

The operator T is not unique (see the remark at the end of this section). 
We define the non-resonant operator Af and the resonant operator 1Z as 

A/"={#o}r, 
n = i-M, 

where the operator 1 is the identity. We notice that the range Rg 1Z of the operator 
1Z is included in Ker{H }. A consequence is that 1ZV is constant under the flow of 
H . 

Let us now assume that Hq is integrable with action-angle variables (A, ip) e 
IxT 1 where IB is an open set of M L and T L is the L-dimensional torus. Thus 
H = H (A) and the Poisson bracket {H,H'} between two elements H and H' of 
A is 

_ dH dIT__dH_ dLT_ 
^ ' ' ~ dA ' ~ Ify ' ~dA' 
The operator {H } acts on V expanded as follows 



as 



V=J2 ^k(A)e ik ^, 
{H }V(A, <p) = J2 ^( A ) ' k ^(A)e ik ^, 



where 

/ * \ dH 

A possible choice of T is 

kez L v ' 

^(A)-k^O 

We notice that this choice of T commutes with {H }. For a given V E A, 1ZV is 
the resonant part of V and AfV is the non-resonant part: 

TIV= J2 Vic(A)e ik ^, (1.4) 

w(A)-k=0 

A/V = J] ^k(A)e ik ^. (1.5) 

oj(A)-k^0 



From these operators defined from the integrable part H , we construct a control 
term for the perturbed Hamiltonian H + V where V £ A, i.e. / is constructed 
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such that H + V + f is canonically conjugate to H + TZV . We have the following 
equation 

e™(H + V + f) = H + KV, (1.6) 

where 

1 - e-i rv } 

f = e-^TZV + A/V - V. (1.7) 

We notice that the operator (1 — e~^ rv ^) / {TV} is defined by the expansion 

l_ e -{rv} °^ (_!)n 

- — = V 1 ' {tvy 

{TV} ^(n + 1)! 1 1 
The control term can be expanded in power series as 

-i)' 



n=l 



We notice that if V is of order e, / is of order e 2 . Therefore the addition of a well 
chosen control term / makes the Hamiltonian canonically conjugate to H + 1ZV. 

If H is non-resonant then with the addition of a control term /, the Hamiltonian 
Hq + V + / is canonically conjugate to the integrable Hamiltonian H + 1ZV since 
TZV is only a function of the actions [see Eq. (II .4ft ] . If H is resonant and 1ZV = 0, 
the controlled Hamiltonian H = H + V + f is conjugate to H . In the case TZV = 0, 
the series ( 11.8ft which gives the expansion of the control term /, can be written as 

oo 
s=2 

where f s is of order e s , where e is the size of the perturbation V, and is given by 
the recursion formula 

fs = -- s {TVJ s . 1 }, (1.10) 

where fi — V. 

Remark : Non-unicity of T- The operator T is not unique. Any other choice 
T' satisfies that the range Rg(r' — T) is included into the kernel Ker({H } 2 ). It is 
straightforward to see that adding a constant to T does not change the expression 
of the control term. However, different control terms are obtained by adding a more 
complex linear operator to T. For instance, we choose 

r = r + kb, 



where B is any linear operator. The operator T' satisfies {Hq} 2 T' = {Hq} since 
{Hq}K = 0. The control term we obtain is given by Eq. (j 1.8 ft where T is replaced 
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by T'. We notice that the new resonant operator defined as TZ' = 1 — {H }T' is 
equal to TZ since {Hq}TZ = 0. The new control term is in general different from /. 
For instance, its leading order is (still of order e 2 ) 

f 2 = f2-\{RBV,{TZ + l)V}. 

This freedom in choosing the operator V can be used to simplify the control term (for 
withdrawing some Fourier coefficients) or more generally, to satisfy some additional 
constraints. 

Remark : Higher- order control terms- It is possible to construct higher order 
control terms, i.e. control terms of order e n with n > 2, and there are many ways 
to do so. In this paragraph, we add e in the control term for bookkeeping purposes. 
For instance, we notice that H^ — H + eV + e 2 f is integrable. The perturbed 
Hamiltonian can be written as 

H = H® - e 2 f. 

For simplicity, we assume a non-resonant condition on V, i.e. TZV = 0. We define a 
new operator F as 

f = e - £{ rv} re£ {ry} ) 
which is e-close to T. It is straightforward to check that T satifies 

{i^ c) } 2 f = {^ c) }, 

since we have {Hq } = e~ e ^ rv ^{Ho}e £ ^ rv ^ from Eq. (11.61) . By applying again 
Eq. (11.61) and replacing V by — e 2 f and H by (and consequently e by — e 2 ), we 
have the existence of a control term e A g which satisfies : 

H + sV + e*g = (^ c) - e 2 TZf) , 

where TZ is the resonant operator associated with f defined as TZ = 1 — {H^}T 
(and is e-close to TZ). It follows from the commutation of {7ZV} and {H^} that 
the controlled Hamiltonian Hq + eV + e A g is integrable. An expansion in e of the 
control term g gives its leading order 

g 4 = -j{T{TV,V},{rV,V}}. 

However, we notice that rjri^, V^} introduces additional small denominators. It is 
not obvious that, given a value of e, the control term e 4 g is smaller than e 2 f (using 
a standard norm of functions). 
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1.3 Reduction of chaotic transport for E x B drift 
in magnetized plasmas 

1.3.1 A paradigmatic model of electric potential 

We consider the following model of electrostatic potential [20] 



V(x, t) = V * sin 
kez 2 



— k ■ x + y? k - u{k)t 



'i.ir 



where tp^ are random phases (uniformly distributed) and Vk decrease as a given 
function of ||k||, in agreement with experimental data [21] , and L is the typical size 
of the elementary cell as represented in Fig. 11.11 In principle one should use for 
oo(k) the dispersion relation for electrostatic drift waves (which are thought to be 
responsible for the observed turbulence) with a frequency broadening for each k in 
order to model the experimentally observed spectrum S*(k, uj). 

Here we consider a quasiperiodic approximation of the turbulent electric potential 
with a finite number K of frequencies. We assume that Uk ^ (otherwise, see remark 
at the end of Sec. 11.3.21) . The phases y?k are chosen at random in order to mimic a 
turbulent field with the reasonable hope that the properties of the realization thus 
obtained are not significantly different from their average. In addition we take for 
||Vk|| a power law in ||k|| to reproduce the spatial spectral characteristics of the 
experimental 5(k), see Ref. [21] • Thus we consider the following explicit form of the 
electric potential: 

K N 

V(x,y,t) = ^2 Y 27 i(n 2 + m 2 ) 3 / 2 sin [ 2n ( nx + my "> + ^ kmn ~ 2nu kt], (1.12) 

k=l m,n=l ^ 1 

where (pkmn are random phases. 

For example, Fig. 11.21 shows contour plots of the potential (I1.12p with only one 
frequency u>± — 1, that is a periodic potential with period 1 in x and y, given by 

a N 1 

V(x,y,t) = — 7-5- sin [2ii{nx + my) + y? nm - 2vrt] , (1.13) 

m.n—l ^ ' 
n 2 -\-m 2 <N 2 

with N = 25. 

Since the model is fluctuating in time, the eddies of Fig. 11.21 are rapidly modified 
in time and where a vortex was initially present, an open line appears, and so on. 
Two particular properties of the model, anisotropy and propagation have been ob- 
served : each image of the potential field shows an elongated structure of the eddies 
and superposing images obtained at different times a slight propagation in the y = x 
direction is found. However, this propagation can easily be proved not to disturb the 
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Figure 1.2: Contour plots of V(x, y, t) given by Eq. (11. 13ft for a — 1, N — 25, t — 
(left panel) and t — 1/2 (right panel), . 

diffusive motion of the guiding centers. The property of propagation can be easily 
understood analytically. In fact, restricting ourselves to the most simplified case of 
an electric potential given only by a dominant mode (n — m — 1) it is immediately 
evident that at any given time the maxima and minima of the sine are located on 
the lines y = — x + constant. As the amplitudes are decreasing functions of n and m, 
this structure is essentially preserved also in the case of many waves. The property 
of anisotropy is an effect of the random phases in producing eddies that are irregular 
in space. 

We notice that there are two typical time scales in the equations of motion: the 
drift characteristic time t^, inversely proportional to the parameter a, and the period 
of oscillation t u of all the waves that enter the potential. The competition between 
these two time scales determines what kind of diffusive behaviour is observed [21 [3] . 
In what follows we consider the case of weak or intermediate chaotic dynamics (co- 
existence of ordered and chaotic trajectories) which corresponds to the quasi-linear 
diffusion regime (see Sec. I1.3.5p . Whereas in the case of fully developed chaos, that 
corresponds to the so-called Bohm diffusion regime, one has to introduce a slightly 
more complicated approach (see remark at the end of Sec. I1.3.2p . 

1.3.2 Control term for a potential which varies rapidly in 
time 

We consider an electric potential of the form V(x, y,t/e), i.e. such that the time 
scale of its variation is of order e. 

We consider a time-dependent Hamiltonian system described by the function 
V(x.,y,t/e), with (x, y) G R 2L canonically conjugate variables and t the time. This 
Hamiltonian system has L + 1 /2 degrees of freedom. 
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We map this Hamiltonian system with L + 1/2 degrees of freedom to an au- 
tonomous Hamiltonian with L + 1 degrees of freedom by extending the phase space 
from (x, y) to (x, y, E, r) where the new dynamical variable r evolves as r(t) = 
t + r(0) and E is its canonically conjugate variable. The autonomous Hamiltonian 
is given by 

H(x,y,T,E) = E + V(x,y,r/e). 
Rescaling r by a canonical change of variable, 

f = r/e, E = eE, 

one obtains 

H(x,y,T,E) = E + eV(x,y,r), (1.14) 

where we have renamed E = E and f = r. In the case of rapidly time-varying 
potentials, Hamiltonian (j!.14p is in the form H = H + eV, that is an integrable 
Hamiltonian H plus a small perturbation eV. 

In our case ifo — E, i.e. independent of x, y, r, therefore we have 

TO = I- 

If we consider a potential V(x, y, r), in the periodic case we can write 

V(x,y,r) = ^^(x,y)e^, 

k 

and the action of T, TZ and A/" on V is given by 

rv = y Yh^ill g*r (L15) 

W = K,(x,y), (1.16) 
J\fV = V(x,y,r) - Vo(x,y). (1.17) 

Otherwise, in the more general case of a non periodic potential, one can write, 
under suitable hypotheses, 

/+oo 
V(K,y,k)e ikT dk, 
•00 

and the action of T, TZ and A/" operators on V is given by 

rv = pv [ +c ° v ^J ,k K ikT dk, (i.i8) 

KV = V(x,y,0), (1.19) 
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and 

A/V = y(x,y,r)-V(x,y,0). (1.20) 
The computation of the control term is now a straightforward application of Eq. (11. 81) . 

Following the previous section, first we map the Hamiltonian system with 1 + 1/2 
degrees of freedom given by Eq. (11.121) to an autonomous Hamiltonian with two 
degrees of freedom. This is obtained by extending the phase space from (x, y) 
to (x, y, E, t) that is considering E the variable conjugate to the new dynamical 
variable r. This autonomous Hamiltonian is 



H(x, y,E,r) = E + V" — — sin[27r(ra + my) + ip kmn - 2iru k T]. (1.21) 



The integrable part of the Hamiltonian from which the operators V, 1Z and M are 
constructed is isochronous 

Ho = E. 

We notice that Hq is resonant (since it does not depend on the action variable x). 
From the action of Y and 1Z operators computed using Eqs. fll.18p - fll.19j> . we obtain 
the control term using Eq. (II. 8ft . For instance, the expression of fi is 



/a(z.y.*)=nZ E 



QfcOfc'( n ' m — nm') 



7T ^— ' u>k{n 2 + m 2 ) 3 / 2 (n' 2 + m' 2 ) 3 / 2 



k'n'm 1 



x {sin[27r((n + n')x + (m + m')y) + <^fcmn + </?*'m'»' ~ 27r(^ fe + a; fe /)t] 

+ sin[27r((n - + (m - m')y) + y?fc m n - y^fc'm'n' - 27r(^ fc - uj k >)t]} (1.22) 

Remark: Similar calculations can be done in the case where there is a zero 
frequency, e.g., u = and u>k ^ for k ^ 0. The first term of the control term is 

f2 = -\{rv,(n + i)v}, 

where 

W = Y] . ° > , sin[27r(nx + my) + <p 0mn ]. 



If we add the exact expression of the control term to H + V, the effect on the 
flow is the confinement of the motion, i.e. the fluctuations of the trajectories of the 
particles, around their initial positions, are uniformly bounded for any time [22] . 
In Sec. 11.3.51 we show that truncations of the exact control term /, like for instance 
f'2 or / 2 + fs, are able to regularize the dynamics and to slow down the diffusion. 

We notice that for the particular model (11.13P the partial control term / 2 is 
independent of time and is given by 
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Figure 1.3: Contour plot of f 2 (left panel) given by Eq. (11.231) and fs at t — for 
a = 1 and iV = 25. 
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Figure 1.4: Contour plot of V + f 2 for a = 0.8 and iV = 25 at t = 1/2. 



_ a 2 ^ n i m 2 - n 2 m 1 

h{X,y,T) - — 2_, ( 7l 2 +m 2 )3 /2( n 2 +m 2 )3 /2 

n 2 ,m 2 

x sin[27r[(ni - n 2 )x + (mi - m 2 )y] + f nimi - fn 2 m 2 }- (1.23) 

Figure [L3] depicts a contour plot of f 2 given by Eq. (11.231) . Figure [L41 depicts a 
contour plot of V + f 2 for a = 0.8 at t — 0. We notice that the controlled potential 
is a small modification of the potential V since this contour plot looks very similar 
to the one decicted in Fig. 11.21 

The computation of the other terms of the series (II. 9p can be done recursively 
by using Eq. (11.101) . Again, for a < 1, V + f 2 + fa is close to V, meaning that the 
control in this regime is based on a small modification of the potential. 
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1.3.3 Properties of the control term 

In this section, we state that for a sufficiently small, the exact control term exists 
and is regular. Then we give estimates of the partial control terms in order to 
compare the relative sizes of the different terms with respect to the perturbation. 
The computations are performed for the model (11.131) but can be easily generalized 
to the more general case (j!.12p . 

Concerning the existence of the control term, we have the following proposition 



Proposition 1 - If the amplitude a of the potential is sufficiently small, there exists 
a control term f given by the series U.9\) such that E + V + f is canonically conjugate 
to E, where V is given by Eq. $1.13\) . 



We have shown [20] that the exact control term exists for a < 7 x 10~ 3 . As 
usual, such estimates are very conservative with respect to realistic values of a. In 
the numerical study, we consider values of a of order 1. 

Concerning the regularity of the control term of the potential (11.131) . we notice 
that each term f s in the series ( 11. 9ft is a trigonometric polynomial with an increasing 
degree with s. The resulting control term is not smooth but its Fourier coefficients 
exhibit the same power law mode dependence as V: 

Proposition 2 - All the Fourier coefficients f^mk °f ^ e functions f s of the se- 
ries U.9\) satisfy: 



I f W I < (i OA) 

for (n, m) ^ (0, 0). Consequently, for a sufficiently small, the Fourier coefficients of 
the control term f given by Eq. ( li.gj) satisfy: 

If I < 

\Jnmk\ _ 



a s C s 



(n 2 + m 2 ) 3 / 2 

for (n, m) ^ (0, 0) and for some constant C^ > 0. 



1.3.4 Magnitude of the control term 

A measure of the relative sizes of the control terms is defined by the electric energy 
density associated with each electric field V, f2 and / 3 . From the potential we get 
the electric field and hence the motion of the particles. We define an average energy 
density S as 

where E(x, y, t) = — W. In terms of the particles, it corresponds to the mean value 
of the kinetic energy (x 2 +y 2 ) (up to a multiplicative constant). For V(x, y, t) given 
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by Eq. flEE]), 



n 2 +m 2 <N 2 

We define the contribution of f 2 and to the energy density by 




87r ' (n 2 + m 2 ) 2 

n,m— 1 ' 



(1.25) 




(1.26) 



For N = 25, these contributions satisfy: 



-| « 0.1 x a 2 . 



It means that the control terms /2 can be considered as small perturbative terms 
with respect to V when a < 1. 

Remark on the number of modes in V: In this section, all the computations have 
been performed with a fixed number of modes N = 25 in the potential V given by 
Eq. (11.131) . The question we address in this remark is how the results are modified 
as we increase N. First we notice that the potential and its electric energy density 
are bounded with N since 



Concerning the partial control term f 2 , we see that it is in general unbounded with 
N. From its explicit form it grows like N log N. Less is known on the control term 
since it is given by a series whose terms are defined by recursion. However, from 
Proposition [21 we can show that the value a of existence of the control term decreases 
like l/(2iVlogiV). This divergence of the control term comes from the fact that the 
Fourier coefficients of the potential V are weakly decreasing with the amplitude of 
the wavenumber. 

Therefore, the exact control term might not exist if we increase N by keeping a con- 
stant. However for practical purposes the Fourier series of the control term can be 
truncated to its first terms (the Fourier modes with highest amplitudes). Further- 
more in the example we consider as well as for any realistic situation the value of N 
is bounded by the resolution of the potential. In the case of electrostatic turbulence 
in plasmas kpi ~ 1 determines an upper bound for k, where k is the transverse 
wave vector related to the indices n, m and pi the ion Larmor radius. The physics 
corresponds to the averaging effect introduced by the Larmor rotation. 





1.3. REDUCTION OF CHAOTIC TRANSPORT FOR E x B DRIFT IN MAGNETIZED PLASM. 



12, 1 1 1 . . 1 12 




Figure 1.5: Poincare sections of a trajectory obtained using a generic initial condition 
for Hamiltonian (I1.13P with a = 0.8, without control term (upper left panel) and 
with control term ( 11. 23ft (upper right panel). 

1.3.5 Efficiency and robustness 

With the aid of numerical simulations (see Refs. [21 [3] for more details on the nu- 
merics), we check the effectiveness of the control term by comparing the dynamics 
of particles obtained from the uncontrolled Hamiltonian and from the same Hamil- 
tonian with the control term f 2 and with a more refined control term f 2 + ^3. 

Diffusion of test particles 

The effect of the control terms can first be seen from a few randomly chosen tra- 
jectories. We have plotted Poincare sections (stroboscopic plots of the trajectories 
of V) on Fig. 11.51 of two trajectories issued from generic initial conditions computed 
without and with the control term f 2 respectively. Similar pictures are obtained for 
many other randomly chosen initial conditions. The stabilizing effect of the control 
term ( 11.231) is illustrated by such trajectories. The motion remains diffusive but the 
extension of the phase space explored by the trajectory is reduced. 

The dynamics is more clearly seen on a Poincare section on the [0, l] 2 torus (i.e. 
by taking x and y modulo 1). Such Poincare sections are depicted in Fig. 11.61 for V 
(upper left panel), V + f 2 (upper right panel) and V + f 2 + f 3 (lower panel). These 
figures shows that the Poincare sections are composed of two types of trajectories : 
ones which are trapped around resonant islands, and diffusive (chaotic) ones which 
lead to global transport properties. We notice that the number of resonant islands 
has been drastically increased by the control term. The mechanism of reduction of 
chaos seems to be that a significant number of periodic orbits has been stabilized 
by the addition of the control term. 

A clear evidence is found for a relevant reduction of the diffusion in presence of 



16 CHAPTER 1. REDUCTION OF THE CHAOTIC TRANSPORT OF IMPURITIES IN TURBULENT MA 




Figure 1.6: Poincare sections on the [0, l] 2 torus for V (upper left panel), V + f'2 
(upper right panel) and V + ji + ^3 (lower panel) with a = 0.6 
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Figure 1.7: On the left panel, mean square displacement (r 2 (t)) versus time t ob- 
tained for Hamiltonian fl 1 . 1 3 1) with three different values of a = 0.7, a = 0.8, a = 0.9. 
In the right panel, diffusion coefficient D versus a obtained for Hamiltonian (11.131) 
(open squares) and Hamiltonian fll . 13j) plus control term (11.231) (full circles). 



the control term (11. 22ft . In order to study the diffusion properties of the system, 
we have considered a set of M. particles (of order 1000) uniformly distributed at 
random in the domain < x, y < 1 at t — 0. We have computed the mean square 
displacement (r 2 (t)) clS du function of time 

1 M 

<^)> = ^Eii x ^)- x '(°)ii 2 

where Xj(i) = (xi(t) , t/i(t)) is the position of the i-th particle at time t obtained 
by integrating Hamilton's equations with initial conditions Xj(0). Figure [L7l (left 
panels) shows (r 2 (t)) for three different values of a. For the range of parameters 
we consider, the behavior of (r 2 (t)) is always found to be linear in time for t large 
enough. The corresponding diffusion coefficient is defined as 

t^oo t 

Figure HTH (right panel) shows the values of D as a function of a with and without 
control term. It clearly shows a significant decrease of the diffusion coefficient when 
the control term is added. As expected, the action of the control term gets weaker 
as a is increased towards the strongly chaotic phase. We notice that the diffusion 
coefficient is plotted on a log-scale. For a = 0.7, the control reduces the diffusion 
coefficient by a factor approximately equal to 10. 
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Figure 1.8: Diffusion coefficient D versus the magnitude of the control term 5 f 2 
where f 2 is given by Eq. (I1.23P for a = 0.7. The horizontal dashed line corresponds 
to the value of D without control term (5 = 0). The dash-dotted line is a piecewise 
linear interpolation. 

Robustness of the control 

In the previous sections, we have seen that a truncation of the series defining the 
control term by considering the first f 2 or the two first terms f 2 + fz in the pertur- 
bation series in the small parameter a, gives a very efficient control on the chaotic 
dynamics of the system. This reflects the robustness of the method. 
In this section we show that it is possible to use other types of approximations of 
the control term and still get an efficient control of the dynamics. 

We check the robustness of the control by increasing or reducing the amplitude of 
the control term [20] . We replace f 2 by 5 ■ f 2 and we vary the parameter 5 away from 
its reference value 5 = 1. Figure [T751 shows that both the increase and the reduction 
of the magnitude of the control term (which is proportional to 5 ■ a 2 ) result in a loss 
of efficiency in reducing the diffusion coefficient. The fact that a larger perturbation 
term - with respect to the computed one - does not work better, also means that 
the control is not a "brute force" effect. 

The interesting result is that one can significantly reduce the amplitude of the control 
(5 < 1) and still get a reduction of the chaotic diffusion. We notice that the average 
energy density e 2 (5) associated with a control term 5 ■ f 2 is equal to e 2 (5) = S 2 e 2 , 
where e 2 is given by Eq. (jl.26p . Therefore, for 5 = 0.5 where the energy neces- 
sary for the control is one fourth of the optimal control, the diffusion coefficient is 
significantly smaller than in the uncontrolled case 5 = (by nearly a factor 3). 
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